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1 Abstract. In this paper we are concerned with a initial boundary-value problem for a 

coupled system of two KdV equations, posed on the positive half line, under the effect of 
a localized damping term. The model arises when modeling the propagation of long waves 
generated by a wave maker in a channel. It is shown that the solutions of the system are 
q . exponential stable and globally well-posed in the weighted space L 2 {e 2bx dx) for b > 0. The 

D 1 stabilization problem is studied using a Lyapunov approach while the well-posedness result 

, is obtained combining fixed point arguments and energy type estimates. 

Oh \ 1. Introduction 

. We consider the initial-boundary value problem for a coupled system of two Korteweg-de 

Vries equation posed on the positive half line under the presence of a localized damping 
represented by a function a, that is, 

u t + u xxx + a 3 v xxx + uu x + aivv x + a 2 (uv) x + a(x)u = 0, x, t £ M+ 
b\v t + v xxx + b 2 a 3 u xxx + vv x + b 2 a 2 uu x + b 2 ai(uv) x + rv x + a(x)v = , . 

u(0,i)=t;(0,t) = 0, t>0 { } 



> 
(N 



(N 



X 



u(x,0) = uq(x),v(x,0) = Vq(x), x € R + , 

where r, ai, 02, 03, 6i, 62 are rea l constants with 61,62 > 0. 

When a = system (jl.ip was derived by Gear and Grimshaw in [9] as a model to describe 
strong interactions of two long internal gravity waves in a stratified fluid, where the two waves 
are assumed to correspond to different modes of the linearized equations of motion. It has 
the structure of a pair of KdV equations with both linear and nonlinear coupling terms and 
has been object of intensive research in recent years (see, for instance, [U [2J El El H21 151 
[lUl [20] ED E21 [23 EDJ). This apparently complicated system also appears as a special case 
of a broad class of nonlinear evolution equations that can also be interpreted as a coupled 
nonlinear version of Korteweg-de Vries generalized equations of the form 

u t + u xxx + f(u,v) x = 
v t + v xxx + g(u, v) x = 0, 

with / and g satisfying f(u,v) = H u (u,v) and g(u,v) = H v (u,v) for a smooth function H. 
The initial boundary value problems for KdV type models arise naturally in modeling small- 
amplitude long waves in a channel with a wavemaker mounted at one end, or in modeling 
coastal zone motions generated by long-crested waves propagating shoreward from deep water 
(see, for instance, Such mathematical formulations have received considerable attention 
in the past, and a satisfactory theory on the global well-posedness of solutions is available in 
the literature. 

Along this work we assume that the real- valued function a = a(x) satisfies the condition 

a G L°°(IR + ) is a positive function and a(x) > ao > a. e. in uj, (1-2) 

where ao is sufficiently large and u contains two sets of the form (0,5) and (/3,oo), for any 
< 5 « 1 and /3 > 0. Observe that in this case the damping term is acting effectively in uj. 
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We also assume that 

r, oi, a 2 , 0.3, b\ and b 2 are real constants with a\ = a 2 , < r, O362 << 1 and b±, 62 > 0. 

Assumption < 0362 < 1 combined with some conservations law satisfied by the solutions al- 
low to obtain a priori estimates leading to the global well-posedness results. The constant r is 
a non-dimensional constant parameter that can be assumed to be very small (see, for instance, 
[3 [29] ). Condition a\ = a 2 is technical and will be used to simplify some computations. 

In order to make more precise the idea we have in mind, let us first consider the total 
energy associated to (jl.ip . in this case 

E(t) = \ ! (b 2 u 2 + b lV 2 )dx. (1.3) 
2 Jr+ 

Under the above conditions, we have 

—E(t) = -[-|«x(0,t) + -v 2 x (0,t) + a 3 b 2 u x (0,t)v x (0,t)} - / a(x)(b 2 u 2 + v 2 )dx 
at 2 2 j K + 

= ~\ \Vhu x (0,t) + fifov x (0,t)) - ~ (1 - a\b 2 ) v 2 x (0,t) - a(x)(b 2 u 2 + v 2 )dx < 0. 

(L4) 

This indicates that terms a(x)u and a(x)v in the equations plays the role of a feedback 
damping mechanism. So, in the light of the computations performed above, we can see that 
(11. ip generates a flow that can be continued indefinitely in the temporal variable, defining 
a solution valid for all t > 0. Moreover, taking (II. 4p into account, the internal stabilization 
problem was addressed by the authors in [22] and a locally exponential decay result was 
derived. In this context, it is important to note that, even when a = or u = 0, the energy is 
decreasing and is strictly so as long as u x (0,t),v x (0,t) 7^ 0. This suggests that any solution 
of the model may converges to zero exponentially in the absence of the internal damping. 
However, this may not always be the case, since the dissipation due to the boundary terms 
does not seem to be strong enough to guarantee the exponential decay of solutions, including 
the case of the scalar KdV equation. In addition, negative result probably holds when the 
support of a is a finite interval. 

Here, two mathematical issues connected to (jl.ip will be addressed; well-posedness and 
large-time behavior of solutions in a weighted space introduced by Kato in [TO]. More pre- 
cisely, for b > we prove that the solutions are globally well-posed and decay exponentially 
to zero in [L 2 (IR + )] 2 , where 

Ll(R + ) = {u : R + -> R; ||ne ba; || L 2 (M+) < 00}. (1.5) 

To our knowledge, both problems have not been addressed in the literature yet and the 
existing results do not give them an immediate answer. As it is known in the theory of 
dispersive wave equations, global well-posedness in any particular space seems to depend 
on the available local theory, on the energy- type inequalities satisfied by the solutions and 
also on the linear theory. Therefore, we first show that the corresponding linear problem 
generates a semigroup of continuous operator in [L 2 (R + )] 2 . In particular, we establish the 
so-called Kato smoothing effect, i. e., the solutions whose initial datum lies in [L 2 (M + )] 2 
not only lies in C([0,T], [L 2 (R + )} 2 ), but also in L 2 (0,T; [i? b 1 (IR + )] 2 ). This property made it 
possible to combine Duhamel formula and a contraction-mapping argument to prove directly 
the local well-posedness result. In order to get the global result we derive some a priori global 
estimate, whence we also deduce the Kato smoothing effect for the nonlinear problem. Those 
a priori estimates were not available. In fact, the only available a priori estimates for the 
initial-boundary value problem (the energy estimate and the kato effect) are not sufficient to 
yield a global well-posedness result. 
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With the well-posedness results in hand, we can investigate the long time behavior of 
solutions. More precisely, we aim to prove that the localized damping mechanism a(x)(u, v) 
is strong enough in order to guarantee the exponential decay of solutions in [L 2 (R + )] 2 , as 
long as the initial data in (jl.ip live inside a ball of radius R (no matter how large is R > 0). 
This is done introducing a convenient Lyapunov function C = C{U{t)) function, which is a 
suitable perturbation of the [L 2 (M + )] 2 -norm of the solution, that satisfies 



where C and u are positive constant with v = v(E(0)). The proof combines energy estimates 
and the result on the exponential decay of E(t) obtained in [22]. It is also worth mentioning 
the work |24| that develops the same analysis for the corresponding scalar KdV equation and 
from which we borrow the ideas involved in our proofs. However, since we are concerned 
with a coupled system of nonlinear equations we need more delicate estimates. Moreover, we 
improve the analysis developed in |24j showing that the inequality above holds without any 
assumption on b. 

Our result on the time decay rate is of local nature in the sense that the exponential decay 
rate is uniform only in bounded sets of initial data. However, the results obtained in this 
paper do not provide any estimate on how the decay rate depends on the radius R of the ball. 
This has been done for nonlinear models, as far as we know, in very few cases and always 
using some structural conditions on the nonlinearity. 

We finish this section mentioning some results obtained for the scalar KdV equation in 
connection with the analysis developed here [S HH [HI HH [13 [23l [251 E3 EZl [28]. In what 
concern system (jl.ip , except in |22| , the problem on internal stabilization was addressed only 
on a bounded domain [31 [161 121] . In fact, most of results present in these works are proved 
following the ideas introduced for the analysis of the corresponding scalar case. 

The analysis described above was organized in two sections. In Section 2 we present the 
global well-posedness result in two parts, the first one is devoted to the linear theory. Section 
3 is devoted the stabilization problem. In both sections we split the results into several steps 
in order to make the reading easier. 



The key of the modern analysis of nonlinear, dispersive evolution equations is the linear 
estimates. Therefore, first we address the well-posedness result for linear problem associated 
to the initial-boundary value problem (jl.ip . 

2.1. The linear system. Fix b > 0. Considerations is first directed to the corresponding 
linear system 



C(U(t)) - C(U ) < c£{U ), for some c> 0. 
Indeed, this inequality gives at once the decay 

£(U(t)) < CC(U )e- ut , 



2. Well-posedness 



u t + u xxx + a 3 v xxx + a{x)u = 



hvt + v xxx + b 2 a 3 u xxx + rv x + a{x)v = 0, x, t G R+, 

u(0,t) = v(0,t) = 0, t > 0, 

u(x,0) = uq(x),v(x,0) = vq(x), X > 0. 



(2.1) 



In order to apply the semigroup theory, we introduce the operator 



A:D(A)c[Ll(R+)] 2 ^[Ll(R+)} 2 



defined by 
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with the domain 

D{A) := S.U = ( M G [^(R + )] 2 : G [7fc(M+)] 2 for 1 < i < 3 and 17(0) = j . 

With the notation introduced above, the following holds: 

Lemma 2.1. The operator A defined above is the infinitesimal generator of a continuous 
semigroup of operators (S(t)) t >o in [L 2 (IR + )] 2 . 

Proof. We first make the change of variables p = e bx u and q = e bx v. Then, if u and v solves 
(|2.ip . p and q solve 



(2.2) 



' Pt + (d x - bfp + a 3 (d x - bfq + a{x)p = 
hq t - rbq + rq x + (d x - b) 3 q + a 3 b 2 (d x - b) 3 p + a(x)q = x,t £ 
p(0,t)=q(0,t)=0, t>0, 

p{x,0) =p (x) = e bx u (x); q(x,0) = q (x) = e bx v (x), x > 0. 
Now we introduce the Hilbert space 

X= [L 2 (IR+)] 2 

endowed with the inner product 



pcf>dx + / qifjdx. 



<l )' \ <■• ))x h 
Next, we consider the operator 

B : D(B) Cl4l 

with the domain 

D{B) = jw = ( M € [# 3 (M+)] 2 : W(0) = 

defined by 

/ \ / — (<9a; — ^) 3 p — 03 (9a: — b) 3 q — a(x)p 

BW = B[ P )=[ r a 3 b 2 3 3 1 a(x) I. (2.3) 

\Qj \-{d x -b)—q-—{d x -b) 6 p-{d x -b) 6 — q--~'' 1 

Under the above conditions, the adjoint of the operator B is the operator B* given by 
/ \ j {d x + bfp + a 3 {d x + 6) 3 g - a(»)p 

where 

B* : D{B*) Cl4l 

and 

D(B*) = \ w = ( P ) G [77 3 (M+)] 2 : W(Q) = d x W(0) = 



Since B is a closed operator (£?** = .B) and the domains D(B) and D(B*) are dense in X, 
the result is proved if the operators B — XI and its adjoint B* — XI are both dissipative for 



some A G M + . Indeed, first pick W 



P 



6 -D(-B) to obtain 



= 7- -fa f PxxxPdx + 3bb 2 [ Pxxpdx - 3b 2 b 2 [ p x pdx + b 3 b 2 [ p 2 dx 
°i I Jr+ Jr+ Jr+ Jr+ 

-a 3 b 2 / q XX xPdx + 3ba 3 b 2 / q xx pdx - 3b 2 a 3 b 2 / q x pdx + b 3 a 3 b 2 / pqdx 

JM.+ JR+ JR+ JR+ 

br f n , 1 



IR+ 

a 3 b 



a(x)p dx 



pq x dx + 
oi Jr+ oi 



bi 



a(x)q dx 



— / (fez - 36fe + 36 2^ - 6 - — / (q xxx - 3bq xx + 3b q x - b q)qdx. 

i Jr+ oi j r + 



(2.5) 



We observe that 



/ 

JR- 



-I'll PxxxPdx = -^pl(Q,t) 



3bb 2 / p xx pdx = —3bb 2 I p x dx 



3b b 2 / p x pdx = 



3ba 3 b 2 / q xx pdx = -3ba 3 b 2 / q x p x dx 



-3b 2 a 3 b 2 / q x pdx = 3b 2 a 3 b 2 I qp x dx 



-a 3 b 2 p xxx qdx = a 3 b 2 q xxx pdx + a 3 b 2 p x (0, t)q x (0, t) 

JR+ JR+ 

3ba 3 b 2 / p xx qdx = -3ba 3 b 2 / p x q x dx. 



Similar computations can be done for the remaining terms. Combining (|2.5|) and the above 
estimates, the following holds 



(BW,W) X 
1 



< 



-3b 



+ Wa 3 b 2 \ 



(b 2 pl + q 2 ) dx - 6ba 3 b 2 / p x q x dx + (b 3 + br) / {b 2 p 2 + q 2 ) dx 

JR+ JR+ 

pqdx -r I pq x dx - \ (b 2 pl(0, t) + ql(0, t)) + a 3 6 2 ^(0, t)p x (0, t) 
H Jr+ 1 



(2.6) 
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In order to conclude this step of the proof we estimate terms on the right hand side of (|2.6 
using Holder inequality: 



66a 3 6 2 / PxQxdx 



26 a 3 b 2 / pqdx 



< 36 I / b 2 p x dx + a^b 2 / Q^dx I 

< (6 3 |a 3 |6i)^ / p 2 dx + (b 3 \a 3 \b 2 ) [ q 2 dx 

b\ Jr+ Jr+ 

< [(&1+&2) |a 3 |6 3 ](^ f p 2 dx+ [ q 2 dx] 

\bi Jr+ Jr+ / 



r / pq x dx 



< -l ^ + \w 2 )vj pdx 



a3Vhq x (0,t)^p x (0,t) < ^q 2 x (0,t) + ^pl(0,t). 

Recalling that a 2 b 2 < 1 and r is sufficiently small (in particular, ^ — 36(1 — 0^62) < 0), from 
the above estimates and (|2.6|) . we get 

(BW, W) x 



1 fr 2 . 



q 2 dx + 



h (b 3 + br)+b 3 \a 3 \(b 1 +b 2 ) + -±- 

lo 2 



a\b 2 ) [ 
Jm 

+ [6 3 + 6r + 6 3 |a 3 |(6i + 6 2 )] I q 2 dx - \ (l - a 2 b 2 ) g 2 (0, t)\ 



62 
61 



p 2 dx 



1 

< — 



(61 + 1) (6 3 + br) + 6 3 |a 3 | (61 + 6 2 ) + £-] ( ^ / p 2 dx + f q' 

*o 2 ] \bi J R + J R + 



dx 



Then, choosing 



(2.7) 



A 



1 



n + 1) (6 3 + br) + 6 3 \a 3 \ (h + 6 2 ) + -f 

ZC»2 



> 



we deduce that 

((B — X)W, W) x <0. 

Hence, B — A is a dissipative in X. Now, picking W G D(B*) and performing as in the 
previous computation, we can deduce that 

(W, = W) x < A (W, W) x , 

which allows us to conclude that B* — A is dissipative. The proof is now complete. □ 

The following existence and uniqueness result for the solution of the system (|2.1|) is a direct 
consequence of Lemma 12.11 and the semigroups theory. 

Theorem 2.2. Let U G [L^(M + )] 2 . Then, there exists a unique weak solution S(t)Uo of 
(jin) suc/i f/wf J7 G C(M + ; [L 2 (R + )] 2 ). If U G -D(A), system /tas a unigue (classical) 

solution U G C(M+; Z>(A)) n C^M+i [.Lg(R + )] 2 ). 



An additional regularity result for the weak solutions of (12. ip is proved in the next Lemma. 
It is essentially the same result proved in Lemma 12.71 

Lemma 2.3. Let U G [L^(M + )] 2 and U the weak solution of (|2.ip . Then, for any T > 0, 
£7 G L 2 (0,T; [i? 1 ^)] 2 ) 

H^llL°°(0,T;[£g] 2 ) + II^IIl 2 (0,T;[L 2 ] 2 ) ^ ^ ll^oll^ap , (2.8) 
where C = C(T) is a positive constant. 



Proof. The proof is obtained following the arguments developed in the proof of Lemma 12.71 
therefore we omit the details. □ 

2.2. The nonlinear system. In order to make the reading easier we split the proof of the 
main result (Theorem 12. 5p into several steps. The first result is taken from [24J. For the sake 
of completeness we repeat it here. 

Lemma 2.4. If U G [F 6 1 (M+)] 2 , then 



i i 



U ^ X [LOO(ffi+)]2 * ^+26111711^^ \\Uf k{m+ , ]2 . (2.9) 



Proof. Let £ G R + and u, u G Hr(R + ). From Cauchy-Schwarz inequality, we get 



'6 1 



u\Oe 2K = J (u 2 e 2bx )jx = (2uu x e 2bx + 2bu 2 e 2bx )dx 
< 2 f / u 2 e 2bx dx] 2 ( f u 2 e 2bx dx] * + 26 / uV^x 



< (2 + 26)||ti|| i a (H+ J|«|| H 



Consequently, 

(n 2 (0 + v 2 (£)) e 2bi < (2 + 26) (||«|| £ 2 (R+) \\u\\ H i {R+) + |M| L 2 (R+) |M|^ 
< (2 + 26) ||^||[ L 2( M+) ]2 \\U\\ [H i bl 



and the result follows. □ 

The main result of this section reads as follows: 

Theorem 2.5. For any Uq G [L 2 (M + )] 2 and T > 0, i/iere exists a unique solution U G 
C([0,T];[L 2 (M+)] 2 )nL 2 (0,T; [^(M+)] 2 ) o/ flU). 

Proof. First, we consider the nonhomogeneous linear problem 

62 ut + h t^xxx + a 3 b 2 v xxx + a{x)u = fi 

hvt + v xxx + b 2 a 3 u xxx + rv x + a(x)v = f 2 , x G R + , t G (0, T), , . 

u (0,t) = «(0,t) =0, ie (0,T), ^- iUJ 
u(x, 0) = uo(x), i>(x, 0) = ^(x), a; G M + . 

For any f 1} f 2 G L 1 ([0, T]; L 2 (K+)) and u , u G L 2 (R+), Lemma ([23]) guarantees that (|2TTDll 
has a unique weak solution, such that u, v G C([0, T]; L 2 (R + )). Moreover, performing as in 
Lemma 12.71 we deduce that 



sup (\\u(t)\\ L2 + \\v(t)\\ Ll ) + ( [ f {u 2 x + v 2 x )e 2bx dxdt 
o<t<r v 6 bJ \h Jr+ 

< C (\\u \\ Ll + \\v \\ q + £ (||/a(i)|L ? + H/2WIL2) df) , 



(2.11) 



for some positive constant C = C(T) nondecreasing in T. 

Now we turn our attention to the local well-posedness of (jl.ip in the space 

F = C([0,T]; [L 2 (M+)] 2 )nL 2 (0,T; [Hl(M+)] 2 ) 
endowed with its natural norm. We introduce the map V defined by 

T(U){t) = S{t)U + [ S{t-s)N(U)(s)ds (2.12) 
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where {S(t)} t >o is the semigroup associated to the linear problem and 

/ -uu x - aivv x - a 2 (uv) x \ 



N(U) 



1 a 2 b 2 aib 2 , , 

i—vv x — uu x — tuv) 

61 bi &i 



(2.13) 



/ 



We shall prove that V has a fixed point in some ball -Br(O) of the space F. Therefore, the 
following result is needed: 

Lemma 2.6. There exists a constant K > 0, such that 

\\T(U 2 ) - T{U±)\\ F < KT\ (\\UxWf + \\U 2 \\ F ) \\U 2 - Ux\\ F , 

V U!,U 2 G F. 



Proof. Let U x 



vi 



U-2 



U 2 
V 2 



G F and T > to be chosen later. Then, using 



(12. lip we obtain a positive constant C > 0, such that 

||r([/ 2 )-r(iq)|| F 



< CL ( \\u 2 U 2x - UlU lx 1 1 £1(0,^1,2) + \\v 2 V 2x - VxVl x \\ L l( 0jT ;L 2 ) + W U2xV2 ~ U ^ V ^ Hl^O.T;^ 



+ \\U2V 2x ~ UlVi x \\ L x^ T . L %} j , 

(2.14) 

1 b 2 

where L = 1 + |ai| + \ a 2 \ H 1 (lad + I a.2 1 ) ■ To estimate the terms on the right hand side 

bi bi 

of the above inequality, we first combine the triangular inequality and Holder inequality to 
obtain 



\u 2 U 2x - U1U1 X \\lI(Q,T;LI) = ll( U 2 - Ul)uix + U 2 (u 2x - Ui x ) 



\LH0,T;I? b ) 



f T u f T „ 

< / \\u 2 — ui\\ LO o ||ui x || L 2 dt + / \\u 2 \\ Loo \\u 2x — ui x \\ L 2 dt 
Jo " b Jo b 

= \\U2 - u i\\l 2 (0,T;L°°) \\ U ^\\l 2 (0,T;L 2 ) + ll U 2x ~ U\x \\ L 2 (0,T;L 2 ) II U 2 1| L 2 (0,T;L°°) 

Analogously, 

\\u 2x V 2 - UlxVl\\ L i^ 0jT . L 2^ = \\u 2x (v 2 - Vi) + Vi (u 2x - Ulx) ||£l( 0) T;£2) 

< \\ u 2x\\l 2 (0,T;L 2 ) 11^2 ~ v l\\ L 2 (0,T;L°°) + ll U 2x ~ u lx\\ L 2 (0,T;L 2 ) IK IIl 2 (0,T;L°°) • 



(2.15) 



(2.16) 



Similar conclusions remains valid for the other terms. Now, we recall that, from Gagliardo- 
Niremberg inequality, we obtain a positive constant Cqn > 0, such that the following in- 
equality holds 



w 



L 2 (0,T;L°°) 



T 



w\\ Loo dt\ < / C GN \\w\\ L 2\\u lx \\ L2 dt 



< C GN [ / \\w\\ L 2 H^lrc Hig dt) <C G N 



\w\\ ,2 dt 



\u± x \\ 2 L 2 dt 



< CgnT* \\M\l°°(p,T;Ll) W U1x Wl 2 (0,T;L 2 ) " 



(2.17) 



Therefore, returning to (|2.14p and taking (|2.15p - (j2.17p into account, we get 
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|r(i7 2 )-r(^ 1 )|| jP 



< CLC GN Ti 



+ "2 



1 1 
\ U iAl 2 (Q,T;L 2 ) \\U2 - «lH2oo(o,T;Xg) H" 2 * ~ Ulx Wh(0,T;L%) 

1 
2 



(0,T;L 2 ) W U2x \\l^(0,T;LI) " U2x Ulx H L 2 (0,T;L 2 ) 



+ vi 



x\\l 2 (0,T;L 2 ) \\ V 2 ~ Vl\\ LOO{0iT . L2) 



\V2x ~ Via 



l L 2 (0,T;L 2 ) 



+ Vo 



L°°(0,T;L 2 ) \\ V2x \\l 2 (0,T-,LI) W V2x Vlx H L 2 (0,T;L 2 ) 



+ [ ll«2x|lLa (0 ,r;^) \\V2 ~ «1 ll L ~(0,T;L?) II ^ ~ ^ H ^(O.T;^) 



+ I ll U lll L -(0,T;L2) INHl2(0,T;L2) H U 2x -Ulx llz^T;^) 

l 1 
+ I ll«2x||x3( ,T;ia) 11^2 - «1 1| ^(o^ig) H U 2z - Ulx||£a ( o,r ; Lg) 



+ «1 



< 



3C LCqnT 



~ (\\ u lx\\ L 2 (0,T;L 2 ) + ll u 2x|| L 2( 0>T;L 2) + IKx || L 2(o iT;L 2) + ||w2x|| L 2( 0iT;L 2) 



+ \\ U i\\l°°(0,T;L 2 ) + |P2||i«>(o,T;L2) + IN \\ L°° (0,T;L 2 ) + ll U 2 Hl<x>(0,T;.L2)J > 

\u 2 x - ^lx|| L 2 (0iT;L 2) + \\v 2x ~ ^Lr|| L 2( 0iT;L 2) + \\u 2 ~ -"l|| L oo( 0jT;L 2) + ||«2 ~ «1 \\ L°° (0,T;L 2 )) 



< KT± {\\U2Wp + ||#i Hi?) ||t^2 — Ux\\p 
The proof of Lemma 12.61 is complete. 



(2.18) 

□ 



The above discussion allows us to conclude the local well-posedness result. Indeed, Lemma 
and Lemma 12.61 show that TU £ F and 



r(U)\\ F <Cl\\Uo\\ [Ll]2 +T* \\U\\ F 



The above estimate guarantee that V maps -Br(O) C F into itself if we choose R = 2C || f^o || [JLg] 2 
and T > sufficiently small. Moreover, for this choice of R and T small enough, Lemma 
2.61 allows us to conclude that T is a contraction in F. Then, by the contraction mapping 
theorem, we obtain the local existence and uniqueness result. The global existence comes 
from the a priori bound of the solutions established in the following Lemmas. 
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Lemma 2.7. Let U E [Lg(R+)] and U the solution of (fTTj) . T/ien, /or any T > 
1 



1 



(&2U (fj^ 1 ) + °i w ( X >T) — b 2 u (x) — biv (x)dx) dx 
+ 2 7 (v^^( .*) + yfi>&2v x (0, *)) dt + ~(l-4b 2 ) v 2 x {0,t)dt (2jg) 
+ [ [ a(x) (b 2 u 2 + u 2 ) dscdi = 



and 



(b 2 u 2 (x,T) + 6iu 2 (x,T) - 6 2 Uo(x) - Mo^)) e 2te dx 

-6oa 3 6 2 / / u x v x e 2hx dxdt + 36 / / (o 2 u 2 + u x) e 2bx dxdt 
Jo Jr+ Jo Jr+ 

-46 3 / T / (b 2 u 2 + v 2 ) e 2bx dxdt - ^ f [ (b 2 u 3 + v 3 ) e 2bx dxdt 
Jo Jr+ 3 Jo J R + 

-8b 3 a 3 b 2 [ [ uve 2bx dxdt-2baib 2 [ [ (uv 2 + vu 2 ) e 2bx dxdt 



(2.20) 



'0 JR+ 

\J (Vhu x (0,t) + Ja^ 2 v x (0,t)\ dt + 1(1 - a 2 o 2 ) ^ v 2 (0,t)dt 

or / T / v 2 e 2bx dxdt + T f a{x) (b 2 u 2 + ?; 2 ) e 26:c dxo!t = . 
Jo ii+ in JR+ 



+ 2 



^4s a consequence, 

II C/ Hl°°(0,T;[L2]2) + ll^|| L 2( 0iT .[ L 2]2) < C ||J7o||[i|]2 , (2.21) 

where C = C(T) is a positive constant. 

Proof. Identity (I2.19P is the classical energy dissipation law taken directly from [22J. It is 
repeated here for the reader's convenience. 

Identity (I2.20p is obtained multiplying the first equation in (jl.ip by b 2 (e 2bx — l) it, the 
second one by {e 2bx — l) v and adding the resulting identities after integrating over M + x 
(0, T). Therefore, the following identities will be needed: 

f {e 2bx -l)uu t dxdt = ^ / {e 2bx -l)u 2 {x,T)dx- h -^ [ (e 2bx - l)u 2 (x)dx (2.22) 
ii Jr+ 2 J R + 2 J R + 

[ (e 2bx - l)u xxx udxdt = 3bb 2 [ [ u 2 x e 2bx dxdt - Ab 3 b 2 [ [ u 2 e 2bx dxdt (2.23) 
Jr+ Jo Jr+ Jo Jr+ 

r ^ [ v x v(e 2bx - l)dxdt = -br C [ v 2 e 2bx dxdt (2.24) 
Jo Jr+ Jo Jr+ 

a2,b 2 I I (e 2bx - 1 ) u xxx vdxdt + a 3 b 2 / / ( e 2bx - 1 ) v xxx udxdt 



JR + ^ ^ JR + (2 ^ 



8o 3 a 3 6 2 f [ uve 2bx dxdt + Qba 3 b 2 / / u x v x e 2bx dxdt. 
Jo Jr+ Jo Jr+ 
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Similar conclusions remains valid for the remaining terms (e 2bx — l)vvt and (e 2bx — l)vv xxx . 
Moreover, since a\ = a 2 , if we put together all the nonlinear terms in (jl.ip . it follows that 



b 2 [' [ u 2 u x (e 2bx - l)dxdt + a x b 2 f 1 
Jo ii+ Jo 



uv)v x (e 2bx - l)dxdt 



+a 1 b 2 [ [ u(uv) x (e 2bx - l)dxdt+ [ [ v 2 v x (e 2bx - l)dxdt 
Jo Jr+ Jo Jr+ 

+aib 2 [ [ (uv)u x (e 2bx - l)dxdt + ai b 2 [ [ v(uv) x (e 2bx - l)dxdt 
Jo Jr+ Jo ii+ 

= b i r / ^ 2hx - ^ dxdt + \ r f ^ x - 

6 Jo ii+ Jo Jr+ 

+aib 2 [ [ (u 2 v) x (e 2bx - l)dxdt + ai b 2 [ [ (v 2 u) x (e 2bx - l)dxdt 
Jo Jr+ x Jo Jr+ x 



[ [ (b 2 u 3 + v 3 )e 2bx dxdt - 2b ai b 2 [ [ (u 2 v + 
Jo ii+ Jo Jr+ 



■lb 

3 



uv 2 )e 2bx dxdt. 



(2.26) 



Then, taking the identities ()2.22[) - (|2.26[) into account, we get 

I f (Jbx _ 1} (6 2U 2 (X)T) _ b2U l( x )) dx + I f (e2to _ i) ( bl v 2 (x,T) - b lV 2 (x)) dx 

* JR+ * Jr+ 

+36 [ [ (b 2 u 2 x + v 2 x )e 2bx dxdt - 46 3 I [ (b 2 u 2 + v 2 )e 2bx dxdt 
Jo Jr+ Jo Jr+ 

+6ba 3 b 2 [ [ u x v x e 2bx dxdt - 8b 3 a 3 b 2 [ [ uve 2bx dxdt 
Jo Jr+ Jo Jr+ 

~— [ T [ (b 2 u 3 + v 3 )e 2bx dxdt-2b ai b 2 [ T [ (u 2 v + uv 2 )e 2bx dxdt 
3 Jo Jr+ Jo Jr+ 



-br [ T [ v 2 e 2bx dxdt+ [ T [ a(x)(b 2 u 2 + v 2 )(e 2bx -l)dxdt = 0. 
Jo Jr+ Jo Jr+ 



(2.27) 



Adding (^T9j) and (12371) . we obtain rf2~20j) . 

Estimate ()2.21j) is proved combining (|2.20p and Gronwall Lemma. Therefore, the next 
steps are devoted to estimate the terms in (I2.20p . From Holder inequality we obtain 



Qba 3 b 2 



< 3b 



T 



u T v T e 2bx dxdt 



o JR+ 
T 



^/m^u x ) haz^/m^v^ e 2bx dxdt 



JR+ 



b 2 u 2 x e 2bx dxdt + 2a 2 3 b 2 / / vle 2bx dxdt] 



o 



(2.28) 



and 



-8b 3 a 3 b 2 



uve 2bx dxdt 



< 4b 3 ( a\b 2 



/ v 2 e 2bx dxdt + 
o Jr+ Jo 



b 2 u 2 e 2bx dxdt 



y 

2.29) 
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In order to estimate the cubic terms we proceed as follows. First, observe that 



/ 

JR- 



u 2 e 2bx dx 



uu T e 2bx dx 



< 



1 



therefore 



/ 

Jr- 



u 2 e 2bx dx < p / 



1 



u 2 e 2bx dx 



.2„2bz, 



u 2 x e 2bx dx 



(2.30) 



Combining Lemma 12.41 and the above inequality we obtain a positive constant c satisfying 



\b 2 u(x)e t ' x \\ L a C(R +- ) < cb 2 \\u x \\ L 2 



(2.31) 



From (|2.3ip . (|2.19p and Holder inequality we deduce that 



2b 

T 



b 2 u 3 e 2bx dxdt < ^ 
o Jr+ 3 



. 2bb 2 „ 



k(x)e fa || L 

i-T 



u 2 e bx dx ) 



T , , f II II f26_ \ (2.32) 

y IFxIIl^ ||n|| L 2 ||n|| L 2 (it < 6 2 y IKIIl 2 ^yC ||u || L 2 ||u|| i2 I (it 



1 1 2 2 
< <™2 11^1^2(0^2) + C1562 ||«||i2( 0>X;i a) 



for any 5 > 0, where c\s = cis(\\Uo\\ 2 L2 p) is a positive constant. Using f|2.32|) the terms that 
are cubic in u, v in (j2.20p can be estimated as follows 

2b ai b 2 f [ (u 2 v + uv 2 ) e 2bx dxdt + ^ ( T ( (b 2 u 3 + v 3 ) e 2bx dxdt 
Jo il+ 3 Jo ./R+ 

(2.33) 



< 25 [b 2 11^1^2(0^2) + |K|| L 2 (0)T;L 2) J + fc 5 [b 2 ||m|| L 2( 0iT;L 2) + ||v|| L 2( 0iT;L 2) J , 

for any 5 > 0, where fe<5 is positive constant that depends on 5 and ||f/o||jx2]2- 

We can now conclude the result. Indeed, identity (|2.20p together with estimates ()2.28|) . 
(^291) . (pT33l) and Holder inequality lead to 

~ [ (b 2 u 2 + M 2 ) e 2bx dx + \ I {b 2 u 2 {x) + hv 2 (x)) e 2bx dx 
2 Jr+ 2 J R + 

< [4b 3 + Ab 3 (a 2 3 b 2 + 1) + br + K s ] [ [ (b 2 u 2 + v 2 ) e 2bx dxdt 

Jo Jr+ 

+25 [ T [ (b 2 u 2 x + v 2 x ) e 2bx dxdt + — [ T [ b 2 u 2 x e 2bx dxdt + 6ba 2 b 2 
Jo Jr+ 2 y Jr+ 



JR+ 



Then, 



vle 2bx dxdt. 

(2.34) 



36 



T 

JR+ 



b 2 u 2 x e 2bx dxdt 



1 [ {b 2 u 2 + b lV 2 )e 2bx dx + 

2 Jr+ 

+ [36 (1 - 2a|& 2 ) - 25] f T [ v 2 x e 2bx dxdt < \ f {b 2 u 2 (x) + hv 2 (x)) e 2bx dx (2.35) 

Jo Jr+ 2 7k+ 

+ [46 3 + 46 3 (a 2 3 b 2 + 1) + br + kg] [ [ (b 2 u 2 + v 2 ) e 2bx dxdt. 

Jo Jr+ 



13 



J 36 36(l-2al& 2 ) I 
Choosing 5 < min < — , > we deduce that 



(b 2 u 2 {x, T) + hv 2 (x, T)) e 2bx dx 



T ( 2 -36) 

<^4^(lKIL ? + IKy+7/ / i b 2U 2 + hv 2 ) e 2bx dxdt, 
2Kb bj J Q J R+ 

where 7 = 7 (||^oll[L 2 ] 2 ) > 0- Then, by Gronwall Lemma we have the following bound for 
the solution 

\\u(T)\\ 2 q + \\v(T)\\ 2 Ll < fa + b 2 ) (\\uo\\ 2 Ll + ll^ollig) e 27T , V T > 0. {2 37) 
Consequently, 

II^IIl°°(o,t ; [l2]2) < C \\Uo\\[ L 2 b] 2 , (2.38) 
where C = C(T) > 0. Now, combining (I2.38P and (12.351) we derive the estimate 

II^IIl 2 (o,T;[l2]2) < C||J7 ||[ L 2]2 , (2.39) 

for some constant C > 0. Finally, from (12T381) and (j!T39j) we obtain (j!T2"T|) . □ 

Lemma 12.71 guarantees that the solution exists globally in time. Uniqueness is standard 
and follows from Gronwall Lemma. Now the proof of Theorem 12.51 is complete. □ 

3. Decay 

The main result of this section is obtained constructing a convenient Lyapunov function 
which is shown to decrease along trajectories. The proof combines the energy identities 
proved in the previous section and the result on the exponential decay for (jl.ip proved by 
the authors in |22j . We observe that, since the total energy E(t) is decreasing, their result 
can be summarized as follows. 

Theorem 3.1. Let a = a(x) be any damping function satisfying (II. 2D . Then, if Uq = 

2 ) e [£ 2 W 

there exists a positive constant C = C(R) satisfying 



€ [L 2 (M)] 2 system (jl.ip is locally uniformly exponential stable, that is, for any R > 0, 



E(0)<C [ \f a(x)(b 2 u 2 + v 2 )dx + l(^u x (0,t) + Ja^b 2 v x (0,t) 
JO l JR+ 1 \ v 



+^(l-o§fe) «*(0,*) 



dt, 



provided E(0) < R. 



Indeed, if the above inequality holds, from (jl.4p we have that E(T) < 7-E(0) with < 7 < 
1, which combined with the semigroup property allow us to derive the exponential decay of 
E(t). 

Now we state our main result. In order to make the reading easier, we split the proof in 
several steps. 

Theorem 3.2. Let U be the solution of (jl.ip given by Theorem \2.5l Then, for any R > 0, 
there exist constants C > and rj > 0, such that 

||f/(t)|| [L 2 ]2 <Ce-^||^o|| [L2]2 , W>0, (3.1) 

whenever 1 1 C/q 1 1 fx, 2 ! 2 — R- 
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Proof. In order to obtain the result, we consider the Lyapunov function 



£(U)W = 11 ( b ^ 2 + M 2 ) e 2bx dx + c b [ (b 2 u 2 + M 2 ) dx, (3.2) 
2 Jr+ Jr+ 



where c b > is constant to be chosen later. Then, to obtain (|3,ip it is sufficient to prove that 

C(U) - C(U ) < -cL(U ), (3.3) 

for some positive constant c. 

We multiply (|2.19|) by 2c b and add the resulting identity to (|2.2U|) to arrive at the relation 

C(U)-C(U ) 

= - Cb J Q (Vhu x {0,t) + J<%b2V x (p,t) \ dt - t%(l - a 2 3 b 2 ) v 2 x {0,t)dt 
- [ [ a(x) (b 2 u 2 + v 2 ) (c b + e 2bx ) dxdt - 6ba 3 b 2 [ [ 

JO Jr+ ^ ' JO Jr~ 



u x v x e x dxdt 



-3b f [ (b 2 u 2 x + v 2 x ) e 2bx dxdt + 4b 3 f T [ (b 2 u 2 + v 2 ) e 2bx dxdt 
Jo Jr+ Jo Jr+ 

I T f ( b2u3 + y3 ) e2bXdxdt + 8b 3 a 3 b 2 [ T [ uve 2bx dxdt 
3 Jo ii+ Jo ii+ 

+2ba x b 2 [ [ (uv 2 + vu 2 ) e 2bx dxdt - \ [ ( y/hu x (0,t) + J ajb 2 v x (0, t)) dt 
Jo Jr+ 2 Jo V v / 

-- (1 - alb 2 ) [ T v 2 x (0, t)dt + br [ T [ v 2 e 2bx dxdt. 
2 Jo Jo ii+ 

Then, from Holder inequality, (f!T28j) . ([2T29]) and (pl33j) we obtain 
C(U)-C(U ) 

<-c fo ^ (y/faMW) + y/a%hv x (0,tfj dt - c b (l - a 2 3 b 2 ) J v 2 x (0,t)dt 

- [ T [ a(x) (b 2 u 2 + v 2 ) (c b + e 2bx ) dxdt + T [ b 2 u 2 x e 2bx dxdt 
Jo Jr+ v ' 2 Jo ii+ 

+66a|6 2 / T / v 2 x e 2bx dxdt + 46 3 a|6 2 / T / v 2 e 2bx dxdt + 46 3 / b 2 u 2 e 2bx dxdt 
Jo Jr+ Jo Jr+ Jo Jr+ 

+k s [ T [ (b 2 u 2 + v 2 ) e 2bx dxdt + 25 [ T [ (b 2 u 2 x + v 2 x ) e 2bx dxdt 
Jo Jr+ Jo ii+ 

+46 3 / T / (b 2 u 2 + u 2 ) e 2fe:E dx^ ~3b [ T [ (b 2 u 2 x + v 2 ) e 2fe:c dxdt 
Jo Jr+ Jo Jr+ 

fy^^(0,t) + ^i^^(0,t)J - i (1 - affc) J v 2 x (0,t)dt 



(3.4) 



+6r / / v 2 e 2bx dxdt. 
Jo Jr+ 



(3.5) 
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Consequently, for any 5 > 0, (|3.5p leads to the inequality 
C(U) - C(U ) 

pT / , \ 2 

< " { c 6 + ' * 



/ T f a(x) (b 2 u 2 + v 2 ) (c,, + e 2bx ^j dxdt - ( y - 25^ / T / b 2 u 2 x e 2bx dxdt 



v 2 x (? hx dxdt 



- (36(1 -2a§6 2 ) -25) / / 

+ [46 3 (1 + a|6 2 ) +br + k s ] [ [ (b 2 u 2 + u 2 ) e 2bx dxdt, 

Jo Jr+ 



(3-6) 

where fc<5 is a positive constant that depends on 5 and ||i7o|| 2 L 2]2- Now, we observe that 



46 3 (l + 03^2) < 66 3 , therefore if we introduce the notation 

p b = [6b 3 + br + kg] , (3.7) 

from (|3.6p we can deduce that 
C(U) - C(U Q ) 

~{ Cb + l)f (^^( ' t ) + '\A#2^( ' t )) dt - ( Cfc + ( l - °3 6 2) J V 2 x {0,t)dt 

[ [ a(x) (b 2 u 2 + v 2 ) (c b + e 2bx ) dxdt + p b [ [ (b 2 u 2 + v 2 ) e 2bx dxdt 
Jo Jr+ ^ ' Jo il+ 



< 



36 
~2 



25) I I b 2 u 2 x e 2bx dxdt - (36 (1 - 2al6 2 ) - 26) I I v 2 x e 2bx dxdt. 
J Jo ii+ " Jo Jr+ 



On the other hand, Theorem 13.11 guarantees the existence of a constant C > satisfying 

nxo pT l-XQ 

(b 2 u 2 + v 2 ) e 2bx dxdt < e 2bxo Pb / / (b 2 u 2 + v 2 ) dxdt 
Jo Jo 

<P&C<~^ (V^u x (0,t) + ^J^ 2 v x (0,t)\ dt + l - (1 - a%) J v 2 x (0,t)dt (3 7) 



+ / / a(x) (b 2 u 2 + v 2 ) dxdt \ . 
Jo it+ 



Moreover, since — a(x) > — ao for x > xo, it follows that 

-/7 ° 

JO JI+ 

*>T /.x /-T r+00 



i(x) (b 2 u 2 + w 2 ) e 2bx dxdt 



nXQ pi p+oo 

a(x)(b 2 u 2 + v 2 )e 2bx dxdt- J J a(x) (b 2 u 2 + v 2 ) e 2bx dxdt ^ 

a{x) (b 2 u 2 + v 2 ) e 2bx dxdt - a / (6 2 u 2 + v 2 ) e 2bx dxdt. 

JO Jx 
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Using ([3]) and (|3|) in ([3]) we obtain the estimate 
C(U)-C(U ) 

1 Pb 



^"^+2 2 



Cj ^ jf (yhUsiO, t) + J$p2V x (0, t)j dt+(l- a 2 3 b 2 ) jf t£(0, 

(c 6 - /)fcC) / / a(x) (b 2 u 2 + v 2 ) dxdt - (a - pb) / (6 2 ^ 2 + v 2 ) e 2bx dxdt 
Jo Jr+ Jo Jx 

[ T r a{x) (b 2 u 2 + v 2 ) e 2bx dxdt - ( ^ - 2^ / T / b 2 u 2 x e 2bx dxdt 
Jo Jo V 2 / Jo ifi+ 

(36 (l - 2a|& 2 ) - 25) / / 



v 2 x e 2bx dxdt, 



from which we conclude that 
C(U)-C(U ) 

fT 



< —A; 

JO JR+ 



jf (^^(0,t) + y^^(0,t)) dt + (1 - a 2 b 2 ) J v 2 x (0,t)dt+ 

a(x) (b 2 u 2 + v 2 ) dxdt + f [ [b 2 u 2 x + v 2 x ) e 2bx dxdt 

Jo Jr+ 



(3.6) 



f 36 36 (l - 2a§6 2 ) 1 

where k is a positive constant. Indeed, first we fix 5 < min < — , > as in the 

proof of (|2.2ip . Then, we consider q, and ao satisfying 

ao > pb and c& > pbC. 
Finally, the proof of (|3.3p is a consequence of © and the following results: 
Lemma 3.3. There exists a constant C > 0, such that 



f £{U){t)dt<C [ [ (b 2 u 2 x + v 2 x )e 2bx dxdt. 
o Jo Jr+ 



(3.6) 



and 



+ (l-a|6 2 ) f vl(0,t)dt+ [ [ (b 2 u 2 x + v 2 x )e 2bx dxdt + [ C(U)(t)dt 
Jo Jo Jr+ Jo 



(3.6) 



Proof. Inequality (j3.3j) is a consequence of f|2.30|) : 



< 



C(U)(t)dt <(c b + -]f [ (b 2 u 2 + b lV 2 ) e 2bx dxdt 

W ( Cb + 2") (1 + 6l) / X+ ^ 2n * + ^ e2tedajd *- 



(3.6) 



To obtain (I3.3j) we proceed exactly as in the proof of (|2.20p . We multiply the first equation in 
(11. ip by (T — t) ue 2bx , the second one by (T — t) ve 2bx and add the resulting identities. After 
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some integrations by parts over M + x (0,T), there appears the relation 

% [ {b 2 u 2 ( X ) + hv 2 (x))e 2bx dx 
^ Jr+ 

= - [ T [ (b 2 u 2 + b lV 2 ) e 2bx dxdt - 36 F [ (b 2 u 2 x + v 2 x ) (T - t) e 2bx dxdt 
2 Jo Jr+ Jo ii+ 

f T f ( ft 2^ 3 + ^ 3 ) (T - t) e 2bx dxdt 
Jo ii+ 

(v^CM + V^^ '*)) (T-t)dt + ^(l-a 2 b 2 ) vl(0,t)(T-t)dt 

-6ba 3 b 2 [ [ u x v x {T - t) e 2bx dxdt + 8b 3 a 3 b 2 [ [ uv{T-t)< 
Jo Jr+ Jo Jr+ 



e 2bx dxdt 



-2b ai b 2 [ [ (uv 2 + vu 2 ) (T - t) e 2bx dxdt 
Jo Jr+ 

+br [ T [ v 2 (T — t) e 2bx dxdt -ft a(x) (b 2 u 2 + v 2 ) (T - t) e 2bx dxdt 
Jo Jr+ Jo Jr+ 

+ f [ a(x) (b 2 u 2 + v 2 ) (T - t) dxdt. 
Jo Jr+ 

(3.6) 

Now, we multiply the identity above by T 1 and observe that < 1. Then, using (jl.2p 
and (|2,33p we can estimate straightforwardly in ([3|) to obtain the inequality 



1 f {b 2 u 2 (x) + b lV 2 (x))e 2bx dx 

* JR+ 



- ( 7 i 7 f + 46 3 |a 3 | 2 6 2 M+l )+/„• + 2 ( 1 + 



So I ( b2u2 + blv2 ) e2bXdxdt 



+ (3b + 25 + 3b\a 3 \ 2 b 2 ) f f (b 2 u 2 x + v 2 x )e 2bx dxdt 

+ (y^ u ^ ^ + \f^v x (0,t)\ di+~(l-a|& 2 ) J v 2 x (0,t)dt. 



From the definition of the function C = C(t), we can conclude that 
1 f {b 2 u 2 (x)+b lV 2 {x))e 2bx dx 

1 JR+ 

~ c '\lo (y^ Ux ^ + ■>/°§k u *(°> t )J dt + (i - a 2 3 b 2 ) v 2 x (o,t)dt 

+ [ [ {b 2 u 2 x + v 2 x )e 2bx dxdt+ [ C(U)(t)dt\ , 
Jo Jr+ Jo J 

where C is a positive constant. Since 

C(U ) < Q + t^j {b 2 u 2 (x) + hv 2 (x)) e 2bx dx, 



(3.6) 



(3.6) 
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the result follows. □ 
Now the proof of Theorem (|3.2p is complete. □ 
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